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Abstract 

The leading and next to leading terms of the average arithmetic area (S(m)} enclosed 
by m — > oo independent closed Brownian planar paths, with a given length t and starting 
£Nj ■ from and ending at the same point, is calculated. The leading term is found to be (S(m)} ~ 

_h | ^ In m and the 0- winding sector arithmetic area inside the m paths is subleading in the 

asymptotic regime. A closed form expression for the algebraic area distribution is also 
obtained and discussed. 
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1 Introduction 
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The question of the area, be it algebraic A or arithmetic S, enclosed by a closed Brownian 
planar path, of a given length t and starting from and ending at a given point, is an already 
old subject. It started in the mid-twentieth century with the well-known Levy's law pQ for the 
probability distribution P{A) of the algebraic area. The question of the probability distribution 
P{S) of the arithmetic area is a much more difficult issue pertaining to its non local nature. 
Interestingly enough, the complete calculation of the first moment of the distribution, (S) = 
irt/5, has been only recently achieved by SLE technics [2]. Since on the other hand the average 
arithmetic area (S n ) = t/{2itn 2 ) of the n ^ 0-winding sectors inside the curve has been known 
for some time [3] thanks to path integral technics, one can readily deduce from these results 
that the average arithmetic area (So) of the 0-winding sectors^ inside the curve is 7rt/30. A 
n-winding sector is defined as a set of points enclosed n times by the path and a 0-winding 
sector is made of points which are either outside the path, or inside it but enclosed an equal 
number of times clockwise and anti-clockwise. Note finally that, having in mind more simple 
winding properties, the asymptotic probability distribution at large time of the angle spanned 
by one path around a given point is also known [3]. 



Clearly, the random variables S n and Sq are such that 5 = J^^L-oo Sn and A = Y^=-oo n ^n- 
They happen to be the basic objects needed to define quantum mechanical models where 
random magnetic impurities are modelised by Aharonov-Bohm vortices [6j. Recently some 
progresses have been made on the geometrical structure of the n-winding sectors thanks to 
numerical simulations of the Hausdorff dimension of their fractal perimeter [7|. 

The question addressed in the present work concerns the generalisation of (S) = irt/5 to 
the arithmetic area (S(m)) spanned by m independent closed paths of a given length t starting 
and returning at the same point. More precisely, one would like to have some information on 
the scaling of (S(m)) when m — > oo. Using again path integral technics in the line of [3], one 
will show that it is possible to compute exactly the leading and next to leading asymptotic 
terms of {S(m) — So(m)) where the contribution of the 0-windings sectors inside the m-paths 
has been substracted precisely for the reason discussed above in the one path case. One will 
find that the leading asymptotic term scales like mm, namely (S(m) — So(m)) ~ y lnm with, 
as already said, no information so far on (So(m)) inside the paths. 
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4 The path integral approach diverges for n — because it cannot distinguish zero- winding sectors inside the 
path from the outside of the path, the latter being of infinite area and, trivially, also a 0-winding sector. 
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To go a little bit further, one might consider a simplification of the problem by looking 
at the arithmetic area of the convex envelop of the paths, a simpler geometrical object than 
their actual fractal envelop. It is known [8J that the area of the convex envelop of one path 
is (S) convex = vrt/2. Recently, the asymptotic behavior for m paths when m — > oo has been 
found [9] to be precisely (S(m)) convex ~ y lnm. In light of the identical scaling obtained 
for (S(m) — So(m)), it means that, in geometrical terms, the m paths tend to fully occupy 
the area available inside their convex envelop when m is large. Clearly, (S(m)) conveic being 
by construction bigger than the actual (S(m)), the latter lies between (S(m) — So(m)) and 
(5(m)) convex . Now, since the leading asymptotic behavior of both the lower and upper bounds 
are found to be equal, in the asymptotic regime one concludes necessarily that (S(m)} ~ ^ mm 
and that (5o(m)) is subleadinj£]. 



2 Algebraic area distribution 

As a warm up, let us consider the generalization of Levy's law to m independent paths. Let 
us first start with one path of length t, starting from and ending to a given point r, so that 
r (0) = r(t) = r. The algebraic area enclosed by this path is A = k. Jq r (r) A df{r) /2 where k 
is the unit vector perpendicular to the plane. The path integral leads to 

,jBAs. _ G B {r,r) 

{ } So(r,r) [) 

where Go(r,r) = l/(2irb) and 

g B (f,r) = / Vr(r)e-ifo ' Md^*-/ = J_ 2 (2) 

Jr(0)=f 27 r*sinhf^J 

In ([1]) the average ( ) has been made over the set C of all paths of length t starting from 
and ending at f. As a result Qb is the Landau propagator of a charged particle in an uniform 
magnetic field B = Bk. Fourier transforming, the probability distribution of A -the Levy's 
law- is 

1 r+00 Bt_ -1 

P(A) = — / 2 e~ iBA dB = (3) 

2^7-00 sinh(f) 2i cosh 2(-4) 

In the case of m paths, one should compute instead of ([3]) 

Bt \ m 



The integral in can be performed by rewriting 1 / sinh(Bf/2) m = 2 m e~ mBt l 2 / '(1 — e~ Bt ) m 
and by expanding, when the integration variable B > 0, the denominator in powers of e *. 
One obtains finally 

Pm(A) = ^ £ ( " + ^ " 1 )( (fe + m/2 | M/t)m+1 + -) (5) 



5 This can be easily understood by noticing that 0- winding sectors tend to disappear when more and more 
paths overlap as m increases, which also implies that n-winding indices tend to increase. 
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Figure 1: The distribution P m (A') for the rescaled algebraic area A' = A/(ty/m). Starting 
from the curve with the greatest maximum m = 1 (Levy's law eq.(J3])), m = 2 (eq. ([?]))■ rn = 3 
(eq. ©), m = oo (eq. ©). 



where the complex conjugate term corresponds to the B < integration (it amounts to set 
A — > —A in the B > integration result). 

Noticing that A is the sum of m independent random variables each satisfying Levy's law, 
one expects A to be gaussian when m becomes large. Indeed, one can observe that, when 
m — > oo, the main contribution to the integral in @ comes from small B values. Thus, 



Bt 
2 



1 



sinM'fW \1 + 



Rescaling the area as A' = A/(ty/m), (|1J) leads to 



Poo (A>) = sT-e-^ 2 (6) 



7T 

In terms of A', we get from (0) or from contour integration, 

p 3 (A') = V^-^ | 3 - 67rv / 3yl'tanh(7rV3yl') - {{^A'f + ^)(1 - 3tanh 2 (7rV3^')) 

2 cosh (7rv3^4') \ 4 / 

(8) 

Figure [T] clearly shows that P m (A') converges quickly to the gaussian P^q^A'^ when vn 
becomes large. In particular A scales like ^/m when m — > oo. The scaling will be different for 
the arithmetic area discussed in the next section. 



3 Winding properties and arithmetic area for m Brownian paths 



3.1 The case of one path: notations and known results. 

The arithmetic area enclosed by a planar Brownian path is closely related to its winding 
properties. Let us again consider a path of length t, starting from and ending at f, and let us 
set 9 to be the angle wounded by the path around a fixed point, say the origin O. Again the 
average {e iae ) is made over C 

i iad\ _ Ga( r i r ) / q \ 



where 

g a (r,r) = Vr(r)e > Jo Jo = — e * JJ V«l J (10) 



By symmetry the average depends only on r with Q a being the propagator of a charged particle 
coupled to a vortex at the origin (the I\ k _ a \s are modified Bessel functions). Obvious symmetry 
and periodicity considerations such as Q a = Q a +i = Qi-a allow to restrict to < a < 1. We 
also set r 2 /t = x (so that 27rrdr = 7rtdx) and (e iae ) = G a (x) with 

+oo 

G a {x) = e- X J2 I\k- a \W (11) 

k=— oo 

so that Gq(x) = 1. 

To integrate over f -while keeping the position of the vortex fixed at the origin- is the same 
as to integrate over the vortex position -while keeping the starting and ending point of the 
path fixed. Therefore it amounts to count the arithmetic areas S n of the n-winding sectors, 
the 0-winding sector included 

roo 

/ irtdxG a (x) = J2(Sny a2 ™ + (S ) (12) 

This integral diverges since, as already stressed in the introduction, S*o is the sum of Sq, the 
arithmetic area of the 0-winding sectors enclosed by the path, and of the area of the outside 
of the path, which is infinite. However, since formally for a = 0, 

roo 

/ irtdxGo(x) = J2(Sn) + {So) (13) 

then r oo 

Z a = Tri / dx (1 - G a {x)) = Y,(S n ) (l - e ,a2ra ) (14) 
is finite. One deduces 

(S n ) = - f 1 Z a e- ia27rn da n ± (15) 
Jo 

and 

(S-S ) = J2{S n ) = J Z a da (16) 
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Using the Laplace transform of the modified Bessel functions, we readily recover 

Z a = nta(l — a) (17) 

(S n ) = ^ n + (18) 

71"/ 

(S-S ) = - (19) 
3.2 The case of m independent paths. 

It is easy to realize that G a (r,r) m provides the appropriate measure for counting the sets of 
m closed paths of length t starting from and ending at r. Following the same line of reasoning 
as in section 13. 1\ we get 

Z a (m) = ntj°°dx(l - G a (x) m ) = £<S„(m)> (l - e ia2 ™) (20) 



and 



(S(m) - S (m)) = T (S n (m)) = [ Z a (m)da (21) 

where S(m), S n (m) and So(m) stand respectively for the total, n- winding sectors and 0- 
winding sectors arithmetic areas enclosed by the m paths. For example a sector of points 
which have been enclosed once by one path in the clockwise direction, twice by another path 
in the anticlockwise direction, and not enclosed by the other m — 2 paths, has winding number 
n = -1 + 2 = 1. 

To find the leading behavior of S(m) in the large m limit, one has to evaluate G a (x) m 
when m — > oo, and so one needs a tractable expression for G a {x). Starting from (jlip and 
using again the Laplace transform of Bessel functions, we get 

dG ^ = A s in(vrQ)e~ 22: (2x) Q - 1 ;7(a - -; 2a; 2x) (22) 

dx yjTT 2 

where U is the degenerate hypergeometric function [TO] 



~—a roo / + \ b— a— 1 



It is easy to verify that G a {x) = G\- a (x), as it should. 
From ([22]) 

G Q (oo) = G a (0) + f°° ^^dx = 1 (24) 
jo dx 

as it should far away from the origin (one has used that when a / 0, G a (0) = 0). More 
precisely, since for z large U(a; b; z) ~ z~ a , one has 

_ . . sin(7ra) _ 2x . 
1 - G a (x) ~ — ^ -e lx when x -)■ oo (25) 

\f2lTX 

On the other hand, for < a < 1/2, 



Changing variable 2x = y In m in ([20 



Z a (m) = ^ Inm f°° dy (l - ( G^V^))™) (27) 
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one sees that in the limit m — > oo 

• when y = 0, that is ylnm = 0, (G a ) m = 

• when y > 0, that is ylnm — )• oo, the asymptotics ([25]) can be used so that 

sin(7ra) j_ 

G a (^r -> fl - ~ e "7^K m (28) 

2 / V ^nylnm J 

It follows that m 

^(^y^)) ^ooe(2/-l) (29) 

where G is the Heaviside function. Finally 

Z a (m) ~ — lnm (30) 

and 

(S(m) — So(m)) ~ — lnm (31) 

As already discussed in the Introduction, it means that since (S(m)) comcx ~ ylnm, at leading 
order 

(S(m)) ~ y lnm (32) 

and, necessarily, (So(m)) is subleading. 

In Figure El numerical simulations (random walks on a square lattice) for (S(m)) corwcx /t, 
(S(m))/t and (S(m) — So(m))/t are displayed. One sees that if (S(m)) convcx /t converges rapidly 
to ^ lnm, this is not the case for (S(m))/t and (S(m) — So(m))/t. It means that subleading 
corrections are needed. They originate from the fact that, when m is large but not infinite, 
(G a ) m deviates from a Heaviside function. One should compute 

1 7T 7T 

-(S(m) - So(m)} = — lnm + — lnm (ci + C2) (33) 



with 



f f y In Tfi N 

C2 = - J da dy \ G a ( — - — ) 



For y > 1, 



so that 



For < y < 1, 



(34) 
(35) 



ylnm A m sin(Tra) 
G a (— — ) ~ m m y (36) 

2 / y 7ry in m 

7T 1 

-lnmci (37) 

^ Vflnm 



, y In m , 
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Figure 2: The average arithmetic area: i) numerical simulations (10000 events) of closed 
random walks of 10 6 steps on a square lattice; m = 4, 8, 16, 1024; +: (S'(m)) convex /t, *: 
{S(m))/t, x: (S(m) — So(m))/t ii) analytical results for (S(m) — So{m))/t; upper continuous 
line: 7iTn(m)/2 (leading order); lower dotted line: eq. (I40p that includes the subleading 
corrections. 



with a = sin(7ra) /\/irylnm ~ sin(7ra) /V 7r In m since when m is large, y is peaked to 1. 
Changing variable to a m}~ y = z and integrating over z, then over a, leads to 



— In m Co ~ In In m (in v4tt — C 

2 2 4 2 V 



+ subleading 



2 " 4 2 V " ' "" / v/7rln(m) 

where C is the Euler constant. Collecting all terms, we finally obtain 

— (S(m) — So(m)) = — lnm — — In In m ^ln\/4vr — C^j + subleading 



(39) 



(40) 



One sees in Figure [2] that the subleading corrections greatly improve the fit: the lower dotted 
line ()40p is indeed not far from the numerical data (the agreement is of course not entirely 
perfect but further subleading corrections seem hard to reach). 



4 Conclusion 



In conclusion one has established that the leading behavior of the average arithmetic area 
(S(m)) enclosed by m — > oo independent closed Brownian planar paths is y lnm. The al- 
gebraic area, on the other hand, scales like t^fm. One should stress that the quite different 
asymptotic behaviors pertain to the essentially different nature of the areas considered: the 
algebraic area is additive, which is not the case of the arithmetic area. On another front, it 
remains a real challenge to get some information on the subleading asymptotic behavior of 
the 0-winding sector area (So(m)). Again path integral technics are not adapted to this case, 
whereas SLE machinery should in principle work. 
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